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ABSTRACT 

In this paper we present the first observations of a large sample of millisecond pul- 
sars at frequencies of 2.7 GHz (Allcm) and 4.9 GHz (A6cm). For almost all sources, 
these represent the first Allcm observations ever. The new measurements more than 
double the number of millisecond pulsars studied at A6cm. Our new flux measurements 



> 

^\j ' extend the known spectra for millisecond pulsars to the highest frequencies to date. 

The coverage of more than a decade of radio spectrum allows us for the first time to 
\ s q ■ search for spectral breaks as so often observed for normal pulsars around 1 GHz. The 

results suggest that, unlike the normal pulsars, millisecond pulsar spectra can be largely 
q^ . described by a single power law. 

We align the observed millisecond pulsar profiles with data from lower frequencies 
Q-f to search for indications of disturbed magnetic fields and attempt to resolve questions 

which were raised in recent literature. Deviations from a dipolar magnetic field structure 
are not evident and absolute timing across the wide frequency range with a single 
dispersion measure is possible. We seem to observe mainly unfilled emission beams, 
which must originate from a very compact region. The existence of non-dipolar field 
K^ , components can therefore not be excluded. 

A compact emission region is also suggested by a remarkably constant profile width 
or component separation over a very wide frequency range. This observed difference 
to the emission properties of normal pulsars is highly significant. For a few sources, 
polarization data at 2.7 and 4.9 GHz could also be obtained which indicate that despite 
the typically larger degree of polarization at lower frequencies, millisecond pulsars are 
weakly polarized or even unpolarized at frequencies above 3 GHz. The simultaneous 
decrease in degree of polarization and the constant profile width thus question proposals 
which link de-polarization and decreasing profile width for normal pulsars to the same 
propagation effect (i.e. birefringence). 
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Comparing the properties of core and conal like profile components to those of normal 
pulsars, we find less significant patterns in their spectral evolution for the population 
of millisecond pulsars. Hence, we suggest that core and conal emission may be created 
by the same emission process. 

Given the small change in profile width, the indicated de-polarization of the radiation 
and the possible simple flux density spectra, MSP emission properties tend to resemble 
those of normal pulsars only shifted towards higher frequencies. 



Subject headings: pulsars: general - radiation mechanisms: non-thermal - stars: evolu- 
tion - polarization 



1. Introduction 

Numerous searches following the discovery of the first and still the fastest rotating neutron star 



PSR Bf937+2f (Backer et al. f982) have revealed about 70 millisecond radio pulsars (e.g. Manch- 
ester 1998). This sample has a characteristic age which is on the average two to three orders of 
magnitude larger than that of normal pulsars, indicating a different evolutionary history. Indeed, 
millisecond pulsars are considered as being spun-up by mass transfer from a binary companion 



( Alpar et al. 1982 ). This evolutionary history distinguishes these recycled pulsars, which we will 
here simply call millisecond pulsars, MSPs, from "normal" pulsars. Chen & Ruderman (1993), for 
instance, suggested that the recycling process has consequences for the topology of the magnetic 
field at the surface of the star which could possibly affect the emission properties of MSPs. 

One of the first investigations of the emission properties of MSPs was carried out by Manch- 
ester (1992). Based upon the small number of MSPs then known, he concluded that they exhibit 
properties similar to normal pulsars. Lorimer et al. (1995) found the spectra of 19 MSPs to be 
slightly steeper than normal pulsars. In a series of papers, we have recently started to study the 
emission characteristics of a large sample of MSPs in a systematic way. Kramer et al. (1998, Paper 
I) found that the spectra of normal and MSPs show great similarities at frequencies up to 1.4 GHz, 
a conclusion which was recently confirmed by Toscano et al. (1998; hereafter TBMS) for a large 
sample of MSPs in the southern hemisphere at frequencies between 400 and 1600 MHz. In Paper I 
we also concluded that MSPs tend to be less luminous compared to normal pulsars. They are also 
found to exhibit narrower emission beams than expected from a period scaling based on data of 
normal pulsars. Xilouris et al. (1998; hereafter Paper II) noted that profiles of MSPs develop much 
less with frequency and often also in an unusual way when compared with a canonical behaviour 



determined from normal pulsars (Rankin 1983, Lyne & Manchester 1988). Studying a very large 



sample of polarization profiles of MSPs for the first time in Paper II, it was also shown that most 
MSP polarization position angle swings appear much flatter than those of normal pulsars. This 
observation was recently confirmed by Sallmen (1998) and Stairs et al. (1999). Their results are 



also in agreement with Paper IPs conclusion that MSP profiles often exhibit a high degree of po- 
larization which remains almost unchanged up to 1.4 GHz (i.e. the highest frequency used in both 
studies). 

Although the data already suggest that the emission mechanism for normal and MSPs are 



essentially identical (see also Jenet et al. 1998), there are prominent differences in their emission 
properties, which could be caused by the different evolutionary history (e.g. the existence of addi- 
tional pulse features in a large fraction of MSPs, see Paper I & II). The most promising approach 
to further investigate this question is a detailed multi-frequency study. In addition to the data 
presented in Paper I and II, recent studies provided new high-quality data for various frequencies 
below 2 GHz, i.e. Sallmen (1998), Stairs et al. (1999) and TBMS. In contrast, observations of MSPs 
at frequencies above 2 GHz are scarce. The examples which can be found in the literature are sum- 
marized in Table |l[ Indeed, prior to the observations presented here, only five MSPs have been 
detected above 3 GHz, viz. PSRs J0437-4715, J1022+1001, J1713+0747, B1855+09 and B1937+21 
(see Tab. |l|). Based on this small sample, Kramer (1998) noted that in various respects the emission 
properties of MSPs up to 5 GHz tend to mimic those of normal pulsars at very high frequencies. 
A study of this apparent trend is essential in understanding the physics of the emission mechanism 
not only in MSPs but in normal pulsars as well. 

The limited number of MSP detections at high radio frequencies prior to this work can be 
mainly attributed to the fact that due to the decrease in flux density at higher frequencies, large 
collecting area and wide observing bandwidths are necessary to obtain sufficient sensitivity to study 
the emission properties in detail. Although larger bandwidths are often available in principle, the 
detection of short-period pulsars over large bandwidths is often hampered by the significant amount 
of dispersive smearing incurred by incoherent detection schemes. Recently, however, a coherent de- 
disperser has been installed at the 100-m radio telescope of the MPIfR in Effelsberg, which provides 
bandwidths of up to 112 MHz. Its combination with a large telescope sensitivity makes this site 
ideally suited for observations of MSPs at high frequencies as shown here. 

In this paper we present the first observations of a large sample of MSPs at frequencies of 2.7 
GHz (Allcm) and 4.9 GHz (A6cm). These are the first Allcm observations for almost all MSPs. 
In this work, we also more than double the number of detections at A6cm. A description of the 
observing system is given in Sect. y. Our new flux measurements presented in Sect. || extend the 
known spectra of MSPs to the highest frequencies studied so far. In Sect. |I| we compare our pulse 
profiles with low frequency measurements in search for similar frequency dependencies as found for 
normal pulsars. After studying the polarization properties for a number of sources up to 5 GHz in 
Sect. ||, we discuss these results and their implications in Sect. ||. 



2. Observations 

All data presented in the following for frequencies above 1 GHz are the result of observations 
carried out between July 1997 to October 1998 with the 100-m Effelsberg radio telescope operated 
by the Max-Planck-Institut fur Radioastronomie, Bonn, Germany. The observing system used for 
1.4 GHz was almost identical to the one described in detail in paper I. However, in contrast to 
Paper I, all profiles presented here were obtained with the Effelsberg-Berkeley Pulsar Processor 
(EBPP), a coherent de-disperser. The EBPP obtains signals of two circular polarizations at an 
intermediate frequency (IF) of 150 MHz provided by the receiving systems and converts these to 
an internal IF of 440 MHz. A bandpass of, at maximum, 112 MHz is split into four portions for 
each of the two circular polarizations, which are mixed down to baseband. Each portion is then 
sub-divided into eight narrow channels via a set of digital filters ( [Backer et al. 1997 ). Their outputs 



are fed into de-disperser boards which in turn perform an on-line coherent de-dispersion. In total 
64 output signals (or 128 if full polarization is desired) are detected and integrated in phase with 
the topocentric pulse period. Sub-integrations of typically 180 s are transferred to the controlling 
computer for further off-line reduction. The total bandwidth of the system is determined by the 
dispersion smearing of the individual source at the observing frequency and can be up to 112 MHz. 
When all four Stokes parameters are requested, the bandwidth of the EBPP is limited to 28 MHz. 

The 2.7 GHz observations employed a cooled HEMT-receiver operating over an 80 MHz band- 
width. The system temperature during the observations was about 40 K on the cold sky with 
a telescope gain of about 1.5 K Jy _1 , comparable to that at 1.4 GHz. At 4.9 GHz we used a 
highly sensitive HEMT receiver with a system temperature of 25 K on the cold sky whilst the 
telescope gain is 1.45 K Jy _1 . This receiver provided either an IF of 750 MHz with a bandwidth 
of 500 MHz or an IF of 150 MHz with a 200 MHz bandwidth. For most of the observations, we 
used the narrow band configuration fed into the EBPP. For some low DM pulsars, we also used 
an incoherent detection system, i.e. the Effelsberg Pulsar Observing System (EPOS), connected to 
a 80-200 MHz multiplying polarimeter in order to obtain an even wider bandwidth than provided 
by the EBPP. A similar set-up was used particularly for polarimetric observations at 2.7 and 4.9 
GHz. For a description of this observing mode we refer to Kijak et al. (1997). The observations 
are summarized in Table ||. 

For comparison with lower frequency data we also present some first results of Effelsberg 
observations made at 820 MHz and 1200 MHz with the EBPP. The HEMT receiver used for these 
observations is tunable between 800 and 1300 MHz and has a system noise temperature of about 
60 K on cold sky, providing a 150 MHz IF connected to the EBPP. Due to the small number of 
observing sessions, reliable flux measurements are not available for these particular observations. 
For all other frequencies, flux density calibrations were performed in the scheme detailed in Paper 
I, i.e. comparing the pulse energy to the output of a switch-able noise diode, which itself was 
calibrated based on observations of standard flux calibrators, such as 3C123, 3C295 and NGC7027 



(|Ott et al. 1994J ). 



The selection of millisecond pulsars to be observed was mainly based on considerations about 
integration time available in a given LST range and the strength of the sources measured at 1.4 
GHz. Guided by our measurements presented in Paper I, we typically selected those sources which 
show peak flux densities of ~20 mJy or larger at 1.4 GHz, i.e. we did not make any selection based 
on spectral index. 



3. Flux densities and spectra at high frequencies 

The spectral indices calculated in Paper I were based on only those high frequency data which 
were available at that time. The extension of the frequency coverage to more than a decade of radio 
spectrum presented in this paper allows us now for the first time to study the spectra in search 
for a steepening in the previously fitted power law spectra as is often observed for normal pulsars 
around 1 GHz (e.g. [Malofeev et al. 1994J ). 



Comparing the results of Paper I to those of TBMS, we find that, with the exceptions of PSRs 
B1620— 26, J1643— 1224 and J2145— 0750, all of the derived spectral indices are consistent within 
the measurement uncertainties. For the three aforementioned pulsars, TBMS's new measurements 
at frequencies below 1 GHz differ from those values in the literature used for Paper I. The low 
frequency data of TBMS represent a very large sample of measurements, and thus they should 
average out scintillation effects more effectively than previous studies might have done. Our flux 
measurements presented in Paper I agree with those of TBMS at the overlapping frequencies, and 
we adopt TBMS's results for the low frequencies. We also include their 1660 MHz flux densities in 
the spectra presented here. For some pulsars, we add flux densities from Stairs et al. (1999) and 
also results of previously unpublished measurements made at a frequency of 820 MHz at the 43-m 
telescope of the National Radio Astronomy Observatory^] in Green Bank. The observing setup and 
calibration procedure of the latter observations, carried out in April 1997 and 1998, are identical to 
those described by Foster et al. (1991). For PSR J2145— 0750 we also include a flux measurement 
made at 102 MHz by Kuzmin & Losovsky (1996). The high frequency flux measurements made 
at Effelsberg typically represent a number of observing sessions of 30-60 min integration time for 
each source. This is sufficient to determine a reliable flux density, since interstellar scintillation 
effects are greatly reduced in amplitude compared to lower frequencies ( Malofeev et al. 1996| ). The 



mean flux densities, S u , are summarized in Table ||. Resulting spectra are presented in Fig. |l| (for 
low frequency data see also Paper I and references therein). Spectral indices were obtained for all 
observed pulsars by a weighted least-squares fit of the power law expression S v oc v a and are quoted 
in Table y. Those sources, which were not detected at 2.7 or 4.9 GHz, but for which interesting 
upper limits on the flux densities can be derived, are collected in Table |I|. 

The first fact which becomes evident from studying Fig. |l] is that most MSP spectra can be 
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essentially fitted by a single power law. This is an interesting result given the variety seen in normal 
pulsar spectra (e.g. Malofeev et al. 1994| ). The fact that a simple power law seems to adequately 



describe MSP spectra well is visible from the difference in spectral indices, Aa, derived for power 
laws separately fitted to data below and above 1.4 GHz (see Table |3|) and from the consistency of 
the spectral indices presented here and in Paper I. Similarly, for pulsars which are common both to 
our and TBMS's study, the spectral indices are in good agreement, even though our uncertainties 
are often smaller due to the larger frequency coverage. The only exception is PSR J1024— 0724, for 
which our observations at 4.9 GHz imply that its spectrum is flatter than found by TBMS (i.e. a 
spectral index of — 1.4±0.1 compared to TBMS's — 1.70±0.12), although it is also well represented 
by a simple power law. 

Given that the spectral indices presented here do not differ significantly from those in Paper 
I for all studied sources, we do not expect a change in any of the conclusions drawn in Paper I. 
Indeed, combining the results of this paper, Paper I and those of TBMS, we still derive a mean 
spectral index for Galactic field MSPs closer than 1.5 kpc of —1.76 ± 0.14. Obviously we have no 
reason to believe that the sample of MSPs studied at high frequencies in this paper comprises a 
special sub-sample of MSPs. In fact, although we did not apply any special criteria to select our 
sources (apart from a modest peak flux density at 1.4 GHz), we detected essentially all MSPs which 
we tried to observe. The derived upper limits for the non-detections listed in Tab. || are consistent 
with the flux densities expected from power law fits made to the lower frequency data. 



4. MSP profiles at multiple frequencies 

In this section we discuss the frequency development of MSP profiles, viz. their relative arrival 
time at various frequencies, the change in profile and component width, the change in component 
separation and also the change in relative component amplitudes. For this purpose we align profiles 
according to their barycentric arrival times at various frequencies using the TEMPO software pack- 



age (Taylor & Weisberg 1989). This time alignment was performed in a similar way as for pulsars 
first described by Foster et al. (1991) and later detailed in Kramer et al. (1997), i.e. constructing 
noise-free templates by describing the profile as the sum of a number of Gaussian components 



(e.g. Kramer et al. 1994). In the ideal case, this procedure is perfectly suited for multi-frequency 
timing measurements. In order to produce a template for an additional observing frequency only 
the relative amplitudes and widths of the components are adjusted while their locations are kept 
fixed. As the chosen fiducial point is therefore the same at all frequencies, a comparison of the 
pulse times-of-arrival (TOAs) at different frequencies is possible in an unbiased way. The TOAs 
are corrected for dispersion delays according to their nominal dispersion measures (see Paper I &: 
II and references therein) and then compared to a timing model for the particular source. During 
this procedure, we usually fit only for the pulse phase, but we also studied the effects of adjust- 
ing the dispersion measure for some cases. Only changes within the nominal uncertainties of the 
known timing solutions were observed, so that dispersion measures were later kept fixed to their 



low frequency values. The results are presented in Figs. [2] - [14], where we add low frequency data 
taken from the European Pulsar Network (EPN) data base which is freely available on the World- 
Wide- Web [] ( [Lorimer et al. 1998 ). References to the authors who generously contributed these 



profiles to the EPN archive, are given in the corresponding figure captions. The profiles of the 
new measurements presented here are also freely available now from the EPN database. We point 
out that for the EPN-profiles timing information was not available. In these cases, alignments 
were done visually by bringing the outer edges of the lower frequency profiles in match with the 
time-stamped data. Before discussing the profile development in general, we make a few notes on 
the individual sources. 

In total we present profiles for fifteen MSPs. For five pulsars this is the first detection around 
5 GHz (i.e. J1012+5307, J1024-0719, J1518+4904, J1643-1224 and J1744-1134) and also for 
nine sources at 2.7 GHz (i.e. J0621+1002, J0751+1808, J1012+5307, J1518+4904, J1640+2224, 
J1643-1224, J1713+0747, J1744-1134 and J2145-0750). Among the eleven MSPs now detected 
at a frequency as high as 5 GHz, PSR J0437— 4715 is the only other source, which is not visible 
from Effelsberg. In order to complete this study, we therefore include data taken from the EPN 
data base provided by Bell et al. (1997) and Manchester & Johnston (1995) (see Fig. |3j) . For all 
sources but PSRs J0751+1807 and J1640+2224 high radio frequency profiles are aligned with data 
obtained at lower frequencies. For the latter two MSPs low frequency data were not available, 
so that we present their 2.7 GHz profiles separately in Fig. [0]. We also include PSR J0218+4232 
here, since our recent measurement at 1.41 GHz represents the highest frequency for which a profile 
is available. Since the number of observations is too small to exclude scintillation effects for this 
source, we do not quote a flux density for PSR J0218+4232, but we align its profile visually with 
data from the EPN archive provided by Stairs et al. (1999). Finally, for PSRs B1620-26 and 
J2051— 0827 we present only spectral information, since the obtained data are sufficient for reliable 
flux measurements but the low signal-to-noise ratio (S/N) of the profiles does not add any additional 
knowledge to their profile development with frequency. All 870, 1220 and 1410-MHz profiles aligned 
with the high frequency data were measured in Effelsberg with the EBPP. The 575-MHz profile 
of PSR J1012+5307 was obtained at Green Bank with the Green-Bank-Berkeley-Pulsar-Processor 
dSallmen 1998|) , which is identical to the EBPP. 



In the following we briefly comment on each individual source observed. For this discussion we 
do not arrange the sources in the usual order of increasing right ascension, but will sort the objects 
into three groups, i.e. isolated, fast spinning and slowly spinning binary MSPs, respectively. This is 
done because in Paper I we have found some indications that profile properties could be correlated 
with the mass of the companion star or, alternatively, with the spin period of the pulsar. Although 
we do not find support for these tentative conclusions of Paper I here, we adopt this scheme for 
the merit of future discussions since it has been suggested (e.g. Chen & Ruderman 1993) that the 
evolutionary history may indeed be reflected at the emission properties. 



The Internet address for this site is http://www.mpifr-bonn.mpg.de/div/pulsar/data/ 



4.1. Isolated MSPs 

In Paper I we confirmed the suggestions made by Bailes et al. (1997) that isolated MSPs tend 
to be less luminous than binary MSPs. Among the eleven sources detected at 4.9 GHz, three MSPs 
are isolated objects, i.e. PSRs J1024-0719 (Fig. 0), J1744-1134 (Fig. 0) and B1937+21 (Fig. ||). 



Pulsar B1937+21 is an extraordinarily luminous source, and PSR J1024— 0719 has a fairly flat 
spectrum as discussed above. By comparison, the spectrum of PSR J1744— 1134 is significantly 
steeper (see Tab. ||), while this source also has a smaller luminosity than the average isolated MSP 
(Paper I). We note that the profile of PSR J1744— 1134 (Fig. [ll]) scarcely changes between 410 
and 4848 MHz, although the presently available S/N prevents a detailed investigation of the high 
frequency data. This is consistent with the pre-cursor, which is detected at 410, 606 and 1408 
MHz, located at a constant separation of about 130° prior to the pulse center. In stark contrast, 
the profile of PSR J1024— 0709 changes dramatically with frequency as visible from the inversion of 
the component amplitude ratio. The central component visible at low frequencies weakens at high 
frequencies, so that only the leading component remains with a weak intensity tail at 4.9 GHz (see 
Fig. 0). 

The profile development of PSR B1937+21 over the whole spectral range is more difficult 
to assess, since at low frequencies the profile is severely scattered by inhomogeneities in the ISM. 
Nevertheless, inspecting the 320 MHz profile presented by Cordes & Stinebring (1984) it seems that 
the trailing component of the main pulse, which is clearly visible at 1.4 GHz, is the dominant feature 
at this low frequency. It is completely disguised by the scattering tail at intermediate frequencies, 
gradually weakens between 1.4, 1.7 and 2.7 GHz and eventually fades towards high frequencies 
(cf. Fig. |l^). At the same time the amplitude ratio of main and interpulse increases gradually 
with increasing frequency. While the first component of the main pulse appears unchanged and 
featureless over the whole spectrum, the second component gradually weakens relative to the main 
pulse (cf. Cordes & Stinebring 1984). This additional component is most likely not associated with 
the observed giant pulses of this pulsar (Jenet et al., in preparation) which appear generally delayed 
relative to the position of the main pulse ([Backer 1995 ; Cognard et al. 1996 ; Sallmen 1998). 



4.2. Fast rotating binary MSPs (P < 10ms) 

The steep spectrum of PSR J0218+4232 combined with its large ratio of dispersion measure 
to period makes it a difficult source to observe at high radio frequencies. Consequently, the 1.4- 
GHz profile shown in Fig. [2] is the highest frequency profile available. Even at 1.4 GHz this pulsar 
still emits over essentially the full pulse period. Although wide profiles are also observed in case of 
normal pulsars (e.g. PSRB1831+04; Lyne & Manchester 1988), the pulse shape of PSR J0218+4232 
and its frequency development are interesting. Indeed, the central component in the main pulse 
feature becomes apparently stronger at higher radio frequencies. At the same time, the weaker 
pulse component preceding the main pulse obviously fades with frequency. Stairs et al. (1999) 



present polarization data of this very weakly polarized pulsar. A determination of the viewing 
geometry is thus difficult, and Stairs et al. (1999) suggest an aligned rotator with a wide beam just 
being grazed. 

The unusually complex profile of PSR J0437— 4715 also develops with frequency (see Fig. || and 
also Fig. 5 of Navarro et al. 1997) in a manner which is different to what is typically observed for 



normal pulsars. In normal pulsars the outer "conal" components become stronger with increasing 
frequency compared to the central "core" component, which eventually even fades away (Rankin 
1983, Lyne & Manchester 1988). In contrast, in the case of PSR J0437— 4715 the amplitude ratio 
of the central core to its outer conal components increases from about 2 to 2.5 at 0.4 GHz to more 
than 13 at 4.6 GHz. Nevertheless, the outer components are still traceable at this frequency, which 
allows us to determine that their separation and the profile width are apparently unaffected by 
a change in frequency. We note however that the position of the core component shifts from a 
slightly off-center position to a more centered (earlier) position at high frequencies. This behaviour 



is comparable to that of the core component of PSR B0329+54 (e.g. Kramer et al. 1997) which 



also shows similarities in the overall pulse shape and polarization characteristics. 

Figure |B| shows a clear profile development for PSR J1012+5307 between 0.6 and 4.9 GHz 
which, however, does not affect the profile width or component separation. The trailing component 
of the main pulse becomes the dominant feature above 1.4 GHz. At the same time, the first 
interpulse feature separated by about 120° from the main pulse brightens to about 80% of the main 
pulse amplitude at 2.7 GHz and is - as far as it can be inferred from the given S/N - as strong as 
the main pulse at 4.9 GHz. This trend has already been noticed in Paper II. 

The small misalignment in Fig. |9| of the profiles of PSR J1643— 1224 between 4.9 GHz and 
the low frequencies is not significant given the relatively low S/N at 4.9 GHz. Overall, the profile 
seems to change very little between low and high frequencies. Indications of a possible scattering 



tail in the 436 MHz profile (Bell et al. 1997) prevent a careful comparison with higher frequency 
data (see also Sect. 4.4). 

The profile of PSR J1713+0747 (Fig. ITH) maintains a highly asymmetric shape throughout the 
whole spectrum. The tail-like feature prominent at low frequencies becomes weaker and almost 
undetectable at high frequencies (see also Kijak et al. 1997). The time alignment suggests a 
simultaneous arrival at all frequencies. 

For PSR B1855+09 a distinct profile evolution has already been demonstrated by Thorsett 
&: Stinebring (1990) and was extensively discussed in Paper II. Here, we only emphasise that the 
separation of the main pulse components shrinks significantly between 430 MHz and 2380 MHz 
(see Fig. 4 of Thorsett fc Stinebring 1990| ). This behaviour is comparable to that for the normal 



pulsars like PSR B1133+16, which also shows the same reverse in amplitude ratio in its outer 



main pulse components as PSR B1855+09 ( phillips fc Wolszczan 19921) . In fact, PSR B1855+09 is 



one of the very rare cases where a MSP shows a frequency development in a manner known from 
normal pulsars (e.g. Rankin 1983 and Lyne &: Manchester 1988). Consistent with that impression 
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is the observation that PSR B1855+09 is also among the very few MSPs, where the position angle 
is not flat but shows an S-like swing (cf. Paper II and references therein). In this context it is 



very interesting to compare the profiles observed at 1.4 GHz and at 4.9 GHz (Kijak et al. 1997) 
as done in Fig. |l^. Relative to the main pulse, the interpulse brightens at high frequencies and 
bifurcates into two distinct components. Therefore, in both main pulse and interpulse the central 
component has a steeper spectrum than the outer components. Applying the knowledge gained 
from normal pulsars (as it is apparently justified for this MSP) this frequency behaviour implies 
that the observed main and interpulse are created at different magnetic poles. In other words, PSR 
B1855+09 is an orthogonal rotator, which is perfectly consistent with the constant separation of 
main and interpulse at all frequencies (see also Foster et al. 1991) and the observed Shapiro delay 
( [Ryba & Taylor 199 1| ; [Raspi, Taylor fc Ryba 1991 ). 



The 2.7-GHz profiles of PSRs J0751+1807 and J1640+2224 shown in Fig. [1J are very similar 
to those presented in Paper I and II for 1.4 GHz. As described in Paper II, the amplitude of the 
first pulse component of PSR J0751+1807 decreases compared to that of the trailing part, while the 
central part seems to disappear completely with increasing frequency. The width of the coherently 
de-dispersed profile presented here is consistent with those of the profiles at 430 MHz and 1410 



MHz ( [Lundgren, Zepka &: Cordes 1995| , Paper I & II) and the coherently de-dispersed one at 1410 



shown by Sallmen (1998). 



4.3. Slow rotating binary MSPs (P > 10ms) 

This group of pulsars is distinct from the previous one by having relatively massive binary com- 
panions, i.e. larger than 0.5 M . All pulsars discussed below are in fact considered to form a single 
family of sources which are spun-up predominantly during a wind-driven accretion phase, i.e. pro- 
ducing a non-monotonic random walk in the angular momentum transferred to the recycled pulsar 
(in contrast to disk accretion with monotonic spin-up; see van den Heuvel 1994) , Arzoumanian, 
pordes fc Wasserman 1999| ). 

The first of this family is PSR J0621+1001 for which we present profiles in Fig. ||. The 
profile remains unchanged between 410 MHz and 2695 MHz and also the component separation 



and component amplitude ratio are remarkably constant (see Sect. 4.4). 



When discussing the profile development for PSR J1022+1001, one has to keep in mind that 
this source is among the sample of MSPs known to exhibit profile instabilities. For a detailed 
discussion we refer the interested reader to Kramer et al. (1999). Pulsar J1022+1001 is another 
example where a MSP exhibits a distinct but unusual profile development with frequency. The 
trailing component, the dominant feature at 370 MHz, is on average weaker than the leading 
component at 430 MHz, but dominates the profiles up 1.4 GHz (cf. Fig. y). Then, towards high 
frequencies, this component weakens again rapidly, being only barely detectable at 4.9 GHz (see also 
Kijak et al. 1997). A similar case is PSR J1518+4904 (Fig. ||) where the intensity of the distinct 
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central component weakens between 370 and 610 MHz, but becomes stronger again at higher 
frequencies, until it is the dominant feature at 4848 MHz. At the same time, the tail-like feature 
disappears similarly as in the PSR J1713+0747 profile. Thus, the leading component appears to 
have a maximum in its flux density spectrum around 400 MHz, which would be consistent with the 
interpretation that the overall spectrum discussed before shows a low frequency turn-over. 

Among the strongest profile developments with frequency is the frequency behaviour of PSR 
J2145— 0750: while the leading component is rather weak at 102 MHz (|Kuzmin fc Losovskii 1996 ), 



it becomes the dominant one at higher frequencies. The trailing pulse peak, strongest at 102 MHz, 
weakens at frequencies of 1.4 GHz and higher, and is hardly visible at 4.8 GHz. The central pulse 
component does not seem to develop at all as it remains the weakest component throughout the 
whole observed spectrum. Unlike the normal pulsars, Fig. [l4| indicates that both the profile width 
and the component separation hardly change between 100 and 4850 MHz. Kuzmin & Losovsky 
(1996) however claim to observe an increase in component separation. The fact that the trailing 
profile feature^] consists of two overlapping components (cf. Paper I and Sallmen 1998) makes a 
clear statement difficult although Sallmen (1998) finds support for a slight increase in component 
separation. An inspection of all coherently de-dispersed profiles available between 410 and 1410 
MHz (Sallmen 1998; Stairs et al. 1999; this work) suggests that neither the component separation 
nor the profile width change significantly (see next section) . Consistently, the prominent pre-cursor 
does not change its position as discussed in Paper II. The clear profile change between 102 and 436 
MHz raises the question as to whether this is common to all MSPs at low frequencies. It would be 
interesting to obtain high quality measurements of a number of MSPs at these frequencies. 



4.4. General frequency behaviour 

From the sample of pulse profiles presented in this paper, two obvious but not necessarily 
expected facts are evident: first, all profiles align very well up to the highest frequencies, and second, 
for some profiles we observe almost no shape changes at all (e.g. PSRs J0621+1002 or J1744— 1134), 
while a few profiles change dramatically which can be most easily described by complete reversals 
in amplitude ratios of profile components (e.g. PSRs J0751+1807 or J2145— 0750). The particular 
examples of PSRs J1518+4904 and J2145— 0750 demonstrate that profile changes might occur at 
frequencies below 400 MHz. However, since such data sets are scarce, we can base any discussion 
only on the profile developments at frequencies of 400 MHz and above. At the high frequency end 
between 2.7 and 4.9 GHz the observed changes are fairly small, although some development can be 
still observed in a few cases (e.g. PSR J1022+1001). In total however, we apparently observe an end 
of profile development around 3 GHz. This frequency seems to be somewhat lower than the typical 



frequency for normal pulsars which lies usually between 5 and 10 GHz (Xilouris et al. 1996, Kramer 



|et al. 1997| ) although it is not inconsistent (cf. Sieber, Reinecke & Wielebinski 1975). Excluding 



6 Note that profiles shown at 102 and 436 MHz have a limiting time resolution of about 10° 
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some possible low frequency profile changes, which may also include scattering effects imparted on 
the pulse profile, the data shown here apparently represent the whole frequency development of 
MSPs in the radio spectrum. 

In this light, it is worth emphasising that there is an obvious lack of a change in component 
separation and MSP profile width. This is a significant difference to the behaviour of normal pulsars. 
In order to demonstrate this behaviour clearly, we followed the example of Foster et al. (1991) by 
measuring the component separation or, if distinct components were not visible, the profile width 
(at a 50% level) as a function of frequency, v, for all MSPs discussed in this work (see Table |6| and 
Fig. [L7]). For measuring profile widths, we followed the example of Paper I and measured a 50% 
intensity level with reference to the peaks of the outermost resolved components. Consequently, 
reference points used in this work and that by Foster et al. (1991) differ for two sources common in 
both samples, PSRs B1620-26 and B1855+09, while for PSR B1937+21 both results are in excellent 
agreement. In addition to the shown profiles we also used data from the literature (e.g. Paper I), if 
appropriate (see Tab. ^), and marked coherently (incoherently) de-dispersed profiles with a circle 
(triangle), respectively. We note that a 50%-width is not an ideal measure in some cases (e.g. PSRs 
J1643— 1224 and J1713+0747 where the existence of additional pulse components at low frequencies 
obviously influence the measurements), but offers the possibility for a direct comparison to results 
published for normal pulsars. For this purpose, the obtained values are modelled to an expression 
W{y) = ao + ai -v 1 by least-squares-fits as done for normal pulsars by Thorsett (1991) and Xilouris 



et al. (1996) (see Fig. 17). Obtained indices, 7, are given in Tab. y. For most sources, the frequency 
dependence is so weak that ao and a\ are highly correlated. We therefore quote a width/separation 
interpolated for 1 GHz instead. We also emphasise that although we measure 50% profile widths 
in some cases and component separations in others, we can readily compare the distribution of 
indices, 7, derived for MSPs and normal pulsars, since both profile width and component separation 



essentially lead to very similar frequency dependencies (cf. Thorsett 1991 and Xilouris et al. 1996). 



The comparison of histograms is made in Fig. 18, which demonstrates very clearly the strikingly 



different behaviour for both types of objects. A Kolmogorov-Smirnov test even yields a probability 
of 10 -8 that both samples are drawn from the same parent distribution. 



5. Polarization 

For the few of the strongest sources at high frequencies or those with a relative low dispersion 
measure-period ratio, we obtained polarization information between 1.7 GHz and 4.9 GHz using 
the EPOS and the EBPP (see Sect. |2|). Observations are summarized in Table || where we quote 
bandwidth, dispersion smearing and the degree of linear and total polarization (i.e. after adding 
linear and circular power in quadrature). The profiles obtained are essentially unpolarized within 
the uncertainties and will be thus shown elsewhere. This result is interesting as many MSPs show 
a modest or even high degree of polarization at lower frequencies (Paper II, Sallmen 1998). For 
comparison we compile data from Paper II, Manchester & Johnston (1995), Navarro et al. (1997), 
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Sallmen (1998), Stairs et al. (1999), Thorsett & Stinebring (1990), and this work, and plot the degree 
of polarization for the linearly polarized intensity as a function of frequency in Fig. [l6|. In Tab. |6] we 
list the de-polarization index, e, which was obtained by fitting the degree of linear polarization as a 
function of frequency to a simple power law, i.e. 11^ oc zA Although the sample is rather small, the 
data suggest that those sources which are highly polarized, de-polarize rapidly, while the weakly 
polarized MSPs de-polarize even further but apparently with a slower rate. We stress that the error 
bars only reflect instrumental uncertainties which at 2.7 and 4.9 GHz nevertheless include estimates 



of a possible de-polarization due to a bandwidth averaging (cf. Melrose &: Macquart 1998J ). This 
is important since measuring the polarization of MSPs apparently involves larger uncertainties 
compared to studies of normal pulsars, since we have discovered in Paper II that the polarization 
characteristics can show significant temporal variations - a result which was later confirmed by 
Sallmen (1998; cf. in particular PSRs J1713+0747 and J2145-0750) and Stairs et al. (1999). As a 
result, the degree of polarization can vary significantly for some sources, e.g. occasionally a profile 
is observed which is much more strongly polarized than the average one. Single pulse studies 
are needed to address this question further, but we believe that an average of many independent 
observation leads to a representative value. 



6. Discussion 

We have demonstrated that the time alignment of MSP profiles presented in our work for 
frequencies between 0.9 and 4.9 GHzy is possible without taking into account timing irregularities 
caused by effects such as aberration, magnetic field sweep-back, magnetic multi-poles or retardation. 
Complementary MSP profiles borrowed from the literature for the sake of completeness of our work 
were aligned by visual inspection with the high frequency profiles. All MSPs have been regularly 
timed at various observatories in the frequency range between 100 and 1400 MHz, which overlaps 
with the frequency coverage in our study. Irregular frequency dependent timing behaviour has not 
been reported so far^] (see e.g. Cordes Sz Stinebring 1984 for PSR B1937+21 as the most strongest 



test). We therefore conclude that, within the resolution and uncertainties of our measurements, we 
do not observe any abnormal frequency dependent timing behaviour over whole the radio spectrum 
observed for MSPs. 



7 0.6 and 4.9 GHz for PSR J1012+5307, respectively. 



8 We note that Stairs et al. (1998) discuss an unusual frequency behaviour of the TO As for B1534+12 which seems 
however to be caused by instrumental effects. 
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6.1. Magnetic field structure 



This absence of any timing irregularities directly reflects on the magnetic field topology in the 
radio emission region, as such could indicate the existence of aberration effects or magnetic multi- 
poles (cf. |Phillips fc Wolszczan 1992 , Kramer et al. 1997 ), Hence, an undisturbed, dipolar structure 



of the magnetic field would be an important result for pulsar emission physics, since the existence of 



non-dipolar field components has been speculated by many authors (e.g. Chen & Ruderman 1993) 
in the past and employed to account for unusual emission phenomena. Among these are (a) the 
smaller than expected emission beams of MSPs in comparison to normal pulsars, (b) the distorted 
polarization position angle curves which often exhibit flat slopes, (c) the abnormal profile frequency 
development and (d) polarization properties (cf. Paper I and II and references therein), or also (e) 
"notches" such as identified in the profile of PSR J0437— 4715 ( fNTavarro fc Manchester 1996| ), 

Usually, field components of higher order than dipolar existing at the surface decay rapidly with 
increasing distance to the neutron star, so that their significance may be low in the radio emission 
region. In addition, there are however two direct contributions to the toroidal component of the 
magnetic field: displacement current (as a rotating dipole in vacuum also has "swept back" field 
lines) and the current of the outsreaming plasma. In MSPs there must be a considerable toroidal 
component present virtually everywhere in the magnetosphere (its relative strength increases as 
{t/Rlc) f° r small r), so that at the light cylinder all the contributions to the magnetic field 
(vacuum dipole poloidal, vacuum dipole toroidal and due to the flowing current) are of the same 
order. 

Given the compactness of MSP magnetospheres these effects could be very well relevant for 
the observed emission properties. For example, the light cylinder radius, Rlc-, which ultimately 
limits the pulsar magnetosphere, is only 70 km for PSR B1937+21. Hence, any non-zero emission 
height may be already at a significant fraction of Rlc- Nevertheless, as we have seen, all MSP 
multi-frequency profiles align perfectly between 0.4 and 5 GHz. This can be interpreted in two 
ways. Either the magnetic field in the whole emission region is purely dipolar, i.e. neither higher 
order multi-poles nor significant toroidal field components are present, or the emission region is so 
compact that we would not notice any deviation from a dipolar structure, i.e. all radio frequencies 
are virtually emitted from the same location. 



6.2. Size of the emission region 

With a typical timing accuracy of about 50/is depending on resolution and S/N, the simultane- 
ous arrival of the multi- frequency profiles sets a limit to the size of the emission region corresponding 
to a light-travel distance of only 2.4 km. Estimating the size of the emission region in this manner 
is rather an oversimplification due to possible non-linear propagation effects in pulsar magneto- 
spheres (e.g. |Barnard &; Arons 1986| , Bcskin, Gurevich fc Istomin 1993| ). Therefore, the size of the 



emission region is also often estimated from the change in profile width (e.g. Cordes 1975 , Phillips 
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Sz Wolszczan 1992, Xilouris et al. 19"96| ) . This approach is based on the observation that the com- 



ponent separation and profile width of normal pulsars typically decreases strongly from low to high 
frequencies (see e.g. [Sieber, Reinecke fc Wielebinski 1975 or Phillips & Wolszczan 1992). This 



behaviour has been explained by the assumption of a frequency dependent location of the emission 
region within the magnetosphere, a so-called radius-to- frequency mapping (RFM, Cordes 1978). 
The location is changing with frequency either in altitude above the surface or in radial distance 
from the magnetic axis. (The latter appears less likely given that profile width and component 
separation change simultaneously.) With the additional assumption that the profile wings reflect 
the extent of the open (dipolar) field line region, the width can be directly converted into the height 
and size of the emission region (cf . Paper I) . 

An alternative explanation for the narrowing of the profile width has been given in terms of a 
propagation effect involving refraction effects and birefringence of the plasma above the polar cap 



(Barnard & Arons 198£, McKinnon 1997, Gallant 1995 , von Hoensbroech, Lesch & Kunzl 1998) 



In this model the emission is created at a single altitude, but not all emission can leave the mag- 
netosphere directly. Instead, one polarization mode first propagates through the plasma before it 
escapes the magnetosphere at a certain radius. The extent of the open field line region at this 
particular, frequency dependent altitude will then be the factor which determines the observed 
profile width. Gallant (1998) gives for the escape radius of the fast O-mode, r esc , an expression 

?~esc „.„/ #12/1472 \ 1/3 -2/3 m 

where R* is the neutron star radius, P the spin period and Byi the surface magnetic field. Interest- 
ingly, unless the multiplication factor between the primary and secondary pair plasma, h = 10 4 /i4, 
and the 7-factor of the secondaries, 7 = 10 2 72, are very low, this escape radius lies outside the 
light cylinder radius for most of the MSPs. If the given approximation is correct, one may expect 
as a consequence to observe only one orthogonal mode of the polarization state. This is generally 
not the case. Even in the case of normal pulsars, the model implies a separation of the orthogo- 
nal modes in pulse phase, but there is clear evidence that the modes occur simultaneously at the 
same phase (e.g. McKinnon fc Stinebring 1998; ). In Paper II this occurrence of orthogonal modes 



was also demonstrated for MSPs, where we observe a discontinuity in the swing of the position 
angle and a drop in the degree of linear polarization for some sources (see also Sallmen 1998 and 
Stairs et al. 1999). Besides, according to the above expression the profile width should exhibit a 
frequency dependence of ex v~ 2 > ?J (modulo some possible geometrical effects). In stark contrast, 
our observations suggest a minimal frequency dependence for profile widths of MSPs. The profile 
width reaches saturation either at very low radio frequencies or does not change with frequency at 
all. 

It is important to emphasise that both models, RFM and birefringence, predict a change in 
profile width with frequency. The difference, however, is that in the RFM model the profile width 
reflects the altitude at which the emission is created, while in the birefringence model it is the 
altitude where the emission escapes the magnetosphere. The birefringence model has certainly the 
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advantage that it can explain the occurrence of orthogonal polarization modes in a natural way. 
Moreover, the RFM model requires the assumption that the overall broad-band pulsar emission 
mechanism exhibits also some narrow-band features in order to produce a distinct mapping of fre- 
quency to altitude, which seems difficult to fulfil simultaneously. Obviously, it appears impossible 
to decide between these models simply on the basis of profile width data. Using additional timing 
data, we can only set an upper limit on the size of the emission region, but the data are in fact con- 
sistent with a single emission altitude (see also Phillips fe Wolszczan 1992 and [Kramer et al. 1997 ). 



Blaskiewicz, Cordes & Wasserman (1991) and von Hoensbroech & Xilouris (1997) have tried to 
determine any dependence of the emission height on frequency by analysing polarization data of 
normal pulsars including relativistic corrections to the rotating- vector- model ( Ftadhakrishnan ~fc| 



Cooke 1969J ). Although this method yields negative emission heights for 5 out of 36 sources (Blask- 



iewicz et al. 1991, von Hoensbroech & Xilouris 1997), and although most of the results are, within 
the (often large) uncertainties, consistent with a constant emission altitude, there is an interesting 
trend visible in the data. According to the interpretation of the Blaskiewicz et al. model, this 
trend suggests a weak but non-zero RFM. However, any evidence has yet to be obtained and a 
single emission altitude is still consistent with the data. The pulsar width data of MSPs presented 
here at least indicate a very weak or even non-existing RFM for these sources, but this may not 
be surprising given the compactness of their magnetospheres — even if RFM is present in normal 
pulsars. 



6.3. De-polarization and profile width 

As mentioned before, the occurrence of orthogonal polarization modes is a natural consequence 
of the birefringence model. An overlap of orthogonal modes leads generally to a de-polarization of 
the average pulse profile. The emission of normal pulsars is in fact well known to exhibit a strong 



de-polarization with frequency (e.g. Manchester, Taylor &: Huguenin 1975) and observations at very 



high radio frequencies (i.e. 32 GHz) suggest that all pulsars become de-polarized at sufficiently high 
frequencies (Xilouris et al. 1996). Among other models (see e.g. Xilouris et al. 1994 for a summary), 



it was therefore suggested that an increased occurrence of orthogonal modes at high frequencies may 
be responsible for the observed de-polarization (|Stincbring et al. 1984). McKinnon (1997) finally 



linked the decreasing profile width and the increasing de-polarization as a function of frequency to 
the same birefringence model. 

Our polarimetric profiles of MSPs discussed in this work suggest that we observe an essentially 
complete de-polarization of MSPs already at frequencies around 5 GHz - at least for the limited 
sample with polarization data up to this frequency. At the same time, even though very little is 
known about single pulse polarimetry of MSPs (see e.g. Sallmen 1998), it is clear that orthogonal 



polarization modes are present in the emission of MSPs (e.g. Thorsett & Stinebring 199C , Paper II, 
Sallmen 1998, Stairs et al. 1999). However, the profiles of MSPs show a strikingly constant width 
across the observed frequencies. In other words, the de-polarization of MSP emission appears to be 
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de-coupled from the effect of profile narrowing. Since MSPs and normal pulsars appear to share the 
same emission mechanism (cf. Sect. |l|), we are therefore led to the conclusion that these phenomena 
are unrelated for normal pulsars as well, at least to zero order. 



6.4. Flux density spectra 

The apparent simplicity of the observed flux density spectra may also point to a small size 
of the emission region, perhaps embedded in a dipolar structure of the magnetic field. One could 
expect that a different field geometry or at least a significant change in emission height also alters 
the efficiency of the emission process, e.g. due to a change in magnetic field strength, curvature 
radius or plasma densities (see e.g. Kuzmin et al. 1986). It is therefore interesting that we very 
often observe a steepening in the flux density spectra of normal pulsars but apparently not for 
MSPs. In fact, at least 55% of all 45 normal pulsars studied by Malofeev et al. (1994) exhibit a 
break in their power law spectra with a mean break frequency of 2.2 ± 0.4 GHz. From the twelve 
MSPs with flux measurements up to 5 GHz, only PSR J0437— 4715 shows a clear steepening in its 
spectrum above 1 GHz, while PSRs J1713+0747 and J2145— 0750 are likely cases. Assuming that 
the fraction of pulsars with double power law spectra is the same for MSPs and normal pulsars, we 
could expect about six MSPs in our sample to exhibit a steepening spectrum. It is important to 
note that the statistical significance for this difference being real may be still low as indicated by 
a x 2 ~t es t which yields only a 60% probability that both samples are different. However, treating 
those spectra which are consistent with a straight power law as such with a break frequency at 
zero Hz, we can compare the spectral break frequency distribution of normal pulsars and MSPs. 
A Kolmogorov-Smirnov test then yields the result, that both samples are drawn from the same 
mother distribution with a probability of only 2.2%. It is also interesting to compare the mean 
spectral index derived for MSPs detected at 4.9 GHz to those of 108 normal pulsars between 1.4 
and 4.9 GHz which were measured with the same observing system by Kijak et al. (1998). The 
MSP spectral index of —1.6 ± 0.1 is somewhat flatter than that of normal pulsars for which the 
mean spectral index is —1.9 ± 0.2 above 1.4 GHz, although we did not apply any spectral criterion 
when selecting our sources. 

An interesting spectral behaviour was already indicated by the early work of Foster et al. (1991) 
who studied the spectra of PSRs B1620-26, B1821-24, B1855+09 and B1937+21. Their analysis 
including also low frequency measurements (cf. Fig. |l|) revealed simple straight power laws from 
very low frequencies up to a few GHz, without indications of any low- frequency turn-over. In fact, 
according to Erickson & Mahoney (1985), there is not any indication of a maximum or turn-over 
in the spectrum of PSR B1937+21 above 10 MHz. The flux measurement at 102 MHz for PSR 
J2145— 0750 by Kuzmin & Losovsky (1996) shown in Fig. |l| establishes this trend for yet another 
source. However, more low frequency flux measurements are of utmost importance to investigate 
this question further. 
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6.5. Profile evolution — Cones & Cores 

We have reviewed in Paper I that the large number of pre- and post-cursors and interpulses in 
MSP profiles is a very significant difference to the emission properties of normal pulsars. In Paper 
II we revealed similarly differences in the profile evolution of MSPs and normal pulsars. These 
findings are confirmed in this work covering a larger frequency range. 

The profile evolution reflects the spectral behaviour of individual profile components. In the 
case of normal pulsars, we observe a very clear dependence of the component spectral index on 
the distance to the magnetic axis ( [Rankin 1983 , Lyne & Manchester 1988 or Kramer et al. 1994). 



Although central components in normal pulsar profiles can be often associated with so-called "core" 
components ( [Rankin 1983 ; Lyne fc Manchester 1988] ), the usually steeper spectral index can also 



be found for some central components which are certainly so-called "conal" components ( |Lyne fc 



Manchester 1988^ Kramer et al. 1994| ). Although we can distinguish core and conal components 



also by their polarisation properties (see [Rankin 1983| ) even in the case of MSPs (see Paper II), 
one should obviously be reluctant to simply identify such obvious differences in emission properties 
directly with different emission processes. Other explanations for the steeper spectrum like a 



geometrical origin are equally likely ( Kramer et al. 1994 , Sieber 1997 ). Whatever is responsible 



for the observed patterns in emission properties of normal pulsars, it is clear from the profile 
developments discussed in Paper II and in this work that the canonical picture of strong central 
components at low and strong outer components at high frequencies does not simply apply to 
MSPs. Although some examples for this trend can be even found among the sample of MSPs 
(e.g. PSRs J0751+1807 and B1913+16), for many sources this is not true and leads either to no 
profile development at all (e.g. PSRs J0621— 1002 and J1744— 1134) or to one which was classified 
as "abnormal" (relative to our expectations originating from studies of normal pulsars) in Paper II 
(e.g. PSRs J0437-4715 or J2145-0750). 

In some cases where we consider the profile development as abnormal, we might have been 
mislead to this classification by the possibility that we observe a profile type which would be 
classified as "partial cone", following the terms of Lyne & Manchester (1988; see their Table 4). In 
such a case, like for the normal pulsar PSR B0355+54, only the leading (or trailing) part of the full 
emission cone is active, so that a core component is located at the trailing (or leading) outer edge 
of the visible profile. As a consequence, this apparent outer (but truly central) component would 
generally show a steeper flux density spectrum. Pulsar J0613— 0200 is almost certainly such a case, 
since the polarization profiles suggest that the trailing component is a core, which also shows by 
far the steepest spectrum (see Sallmen 1998 and Stairs et al. 1999). Therefore, the study of the 
profile development of MSPs shows that unfilled emission beams are certainly present. This fact 
may help to solve the problems of the undersized emission beams discussed in detail in Paper I 
and the overall flat position angles which were first noted in Paper II (see Sallmen 1998). However, 
we are still left with the fact, that typical patterns in profile development are hardly observed for 
MSPs. 
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The lack of clear trends visible in the profile evolution of MSPs (with frequency) is a good 
indication that either the emission structures are not as regular as for normal pulsars (e.g. due 
to unfilled emission beams), or that the formation of recognisable patterns is only possible if the 
emission region provides enough space. Again, the compactness of the magnetosphere might thus 
prevent the formation of profile components with distinct emission properties. The presence of 
disturbing effects like swept-back magnetic fields or magnetic multi-poles is still possible but less 
likely given the good alignment of the pulse profiles across the frequencies. Under the presumption 
that the emission mechanism of MSPs and normal pulsars is the same, as supported by the results 
of Paper I h II and also Jenet et al. (1998), this study of MSPs suggests that the different char- 
acteristics of core and conal components found for normal pulsars are in fact due to differences in 
emission heights or size of the emission region or due to geometrical or propagation effects in the 
magnetosphere, but not due to a fundamentally different emission process. This conclusion is in 
agreement with that of Lyne & Manchester (1988) which recently found new support by the work of 
Manchester et al. (1998), who presented a large number of polarization profiles of normal pulsars, 
where "typical" core polarization properties can be sometimes observed in cone components. 



7. Conclusions 

Considering the frequency independence of the profile width and the observed de-polarization 
of the radiation, MSP emission properties tend to resemble those of normal pulsars only shifted 
towards higher frequencies. This trend is also consistent with the early termination of the profile 
development in frequency. All these characteristics may be understood in the context of a very 
compact magnetosphere existing for MSPs which possibly contains only a extremely small emission 
region. Most likely, these regions are too small to allow for timing irregularities or a significant 
change in profile shape or width or perhaps flux density spectrum. 

While the profile width and component separation appear remarkably constant for all frequen- 
cies, we cannot ultimately decide in favour or against a RFM, which would be certainly weak for 
MSPs. However, polarization observations for a few sources up to 5 GHz suggest that the decrease 
in profile width and de-polarization of pulsar emission are not directly related. 

We find weak indications that the spectra of MSPs may not show the spectral break as often 
observed for normal pulsars around a few GHz. At the same time, we cannot observe a clear 
dependence of spectral index of individual pulse components on distance to the profile midpoint 
(and presumably magnetic axis). Nevertheless, features like partial cones as known from normal 
pulsars seem also to be present for some MSPs. A combination of all available MSP data with 
the previously shown fact that the radiation of MSPs and normal pulsars is obviously created 
by the same emission mechanism, leads us to conclude that there is a no direct evidence for a 
fundamentally different emission process of so-called core and conal components. 

Finally, this study demonstrates that for a number of MSPs regular observations at frequencies 
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above 2 GHz are easily possible with a wide bandwidth system. Besides studying emission proper- 
ties, such measurements can provide a useful tool for high precision timing by a determination of 
disturbing effects due to the interstellar weather ( packer fc Wong 1998 ). 



During this work we made extensive use of the European Pulsar Network data archive. We are 
thus indebted to all authors who generously contributed their measurements to the database. MK 
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from stimulating discussions with J. Arons, J. Hibschman, A. Melatos, S. Sallmen, A. Sorrier, and 
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Fig. 1. — Flux density spectra for 20 millisecond pulsars discussed in this work. For references see 
text. 

Fig. 2.— Profiles of PSR J02 18+4232 at different frequencies. The 410-MHz and 610-MHz profiles 
were taken from Stairs et al. (1999). Profiles are aligned by eye since timing information was not 
available for the low frequency data. 

Fig. 3. — Profile of PSR J0437— 4715 at different frequencies. The alignment of the profiles, which 
were provided to the EPN-data archive by Manchester & Johnston (1995) and Bell et al. (1995), 
was made by eye following the measurements by Navarro et al. (1997). 

Fig. 4.— Profiles of PSR J0621+1002 at different frequencies. The 410-MHz profile was taken 
from Stairs et al. (1999), the 575-MHz profile from Camilo et al. (1996). Profiles are aligned by 
eye since timing information was not available for the low frequency data. 

Fig. 5. — Profiles of PSR J1012+5307 at different frequencies. The 575-MHz profile was taken 
from Sallmen (1998). Profiles are aligned by a timing model, confirming the weak detection at 4.9 
GHz. 

Fig. 6. — Profile PSR J1022+1001 at different frequencies. The 370-MHz profile was taken from 
by Sayer et al. (1997), those at 410 and 610 MHz from Stairs et al. (1999). Apart from these profiles 
the data are aligned by a timing model. 

Fig. 7. — Profiles of PSR J1024— 0709 at different frequencies. The 436-MHz profile was taken 
from the EPN-database and was provided by Bailes et al. (1997). Apart from this profile the data 
are aligned by a timing model. 

Fig. 8.— Profiles of PSR J1518+4904 at different frequencies. The profiles at 370, 575 and 800 
MHz were taken from Nice et al. (1996), the profile at 610 MHz from Stairs et al. (1999). Apart 
from these profiles the data are aligned by a timing model. 

Fig. 9. — Profiles of PSR J1643— 1224 at different frequencies. The profile at 610 MHz was taken 
from Stairs et al. (1999). Apart from this profile the data are aligned by a timing model. 

Fig. 10.— Profiles of PSR J1713+0747 at different frequencies. The 430-MHz and 2305-MHz 
profiles were taken from Foster et al. (1993), the 610-MHz profile from Stairs et al. (1999). Apart 
from these profiles all data are aligned by a timing model. 

Fig. 11. — Profiles PSR J1744— 1134 at different frequencies. The 410-MHz profile was taken from 
Bailes et al. (1997), the 610-MHz profile from Stairs et al. (1999). Apart from these two profiles 
all data are aligned by a timing model. 

Fig. 12. — Comparison of the profiles of B1855+09 at 1.41 and 4.85 GHz. The latter was taken 
from Kijak et al. (1997). Both profiles were aligned by a timing model. 
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Fig. 13. — Profiles of PSR B1937+21 at different frequencies. All profiles are aligned by a timing 
model. 

Fig. 14.— Profiles of PSR J2145-0750 at different frequencies. The profile at 102 MHz was 
provided by Kuzmin & Losovskii (1996), the profiles at 410 and 610 MHz by Stairs et al. (1999). 
The 4850-MHz profile was taken from Kijak et al. (1997). All profiles above 610 MHz are aligned 
by a timing model. 

Fig. 15.— Pulse profiles of PSRs J0751+1807 and J1640+2224 at 2.695 GHz. 

Fig. 16. — Degree of polarization for the linearly polarized intensity. Data were compiled from 
Paper II, Manchester k Johnston (1995), Navarro et al. (1997), Sallmen (1998), Stairs et al. (1999), 
Thorsett & Stinebring (1990) and this work. 

Fig. 17.— Profile width measured at a 50% intensity level (filled symbols) and component sep- 
aration (open symbols) for the sources discussed in this work. Circles mark measurements where 
coherently de-dispersed profiles were used, triangles mark values based on incoherently de-dispersed 
data. References for the used profiles are given in the text, figure captions and Paper I and II. 

Fig. 18. — Comparison of index of profile narrowing (unfilled area) and component separation 
change (shaded area) for normal and millisecond pulsars. The values for normal pulsars were taken 
from Thorsett (1991) and Xilouris et al. (1996). 
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Table 1. Previous observations of millisecond pulsars above 2 GHz. 



PSR 



Period 
(ms) 



Frequency 
(GHz) 



Ref. 



J0437-4715 


5.757 


2.3 


1 






4.6 


2 


J0613-0200 


3.062 


2.3 


1 


J1022+1001 


16.453 


2.3 


3 






4.9 


4 


J1045-4509 


7.474 


2.3 


1 


J1643-1224 


4.621 


2.3 


1 


J1713+0747 


4.570 


2.3 


5 






4.9 


4 


B1855+09 


5.362 


2.3 


6 






3.0 


7 






4.9 


4 


B1937+21 


1.558 


2.3 


7 






3.0 


8 






4.8 


9 






8.5 


9 


J2145-0750 


16.052 


4.9 


4 



References. 



1-Bcllctal. 1997,2-Manch- 



cstcr k Johnston 1995| , 3 - pamilo 1994 4 
Kijak et al. 1997, 5 - Foster, Wolszczan & 



Camilo 199S, 6 



Thorsett fe Stincbring 1990 



7 



Stincbring 1983, 8 - Foster, Fairhead fc 



Backer 1991, 9- Ryba 1991 
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Table 2. Observations of millisecond pulsars presented in this work. 



PSR 


Period 


Frequency 


Backend 


Bandwidth 


Resolution 51 




(ms) 


(GHz) 




(MHz) 


( M s) 


J0218+4232 


2.323 


1.4 


EBPP b 


45 


143 


J0621+1002 


28.853 


1.4 


EBPP 


56 


113 






2.7 


EBPP 


80 d 


197 


J0751+1807 


3.479 


2.7 


EBPP 


80 


31.9 


J1012+5307 


5.256 


0.9 


EBPP 


56 


36.0 






1.2 


EBPP 


56 


36.0 






1.4 


EBPP 


90 


15.4 






2.7 


EBPP 


80 


46.0 






4.8 


EBPP 


112 


77.0 


J1022+1001 


16.453 


1.4 


EBPP 


28 


32.1 






2.7 


EBPP 


80 


64.3 






4.8 


EBPP 


112 


96.4 


J1024-0719 


5.162 


1.4 


EBPP 


45 


55.5 






4.8 


EBPP 


112 


40.3 


J1518+4904 


40.935 


1.4 


EBPP 


28 


80.0 






2.7 


EBPP 


80 


160 






4.8 


EBPP 


112 


160 


J1640+2224 


3.163 


2.7 


EBPP 


80 


27.3 






4.8 


EBPP 


112 


38.6 


J1643-1224 


4.621 


1.4 


EBPP 


23 


18.0 






2.7 


EBPP 


90 


22.6 






4.8 


EBPP 


112 


27.1 


J1713+0747 


4.570 


1.4 


EBPP 


56 


8.9 






1.7 


EBPP 


45 


8.9 






2.7 


EBPP 


80 


17.9 






4.9 


EPOS c 


80 


93.0 


J1744-1134 


4.075 


1.4 


EBPP 


112 


15.9 






2.7 


EBPP 


80 


31.8 






4.8 


EBPP 


112 


31.8 


B1855+09 


5.362 


1.4 


EBPP 


28 


10.5 






4.9 


EPOS 


80 


77.4 


B1937+21 


1.558 


0.9 


EBPP 


28 


1.5 






1.4 


EBPP 


28 


1.5 






1.7 


EBPP 


28 


3.0 






2.7 


EBPP 


80 


6.1 






4.8 


EBPP 


112 


12.1 


J2145-0750 


16.052 


1.4 


EBPP 


90 


31.4 






2.7 


EBPP 


112 


62.7 






4.9 


EPOS 


500 


329.1 



a final resolution after data reduction 
b coherently de-dispersed 
incoherently de-dispersed 
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Table 3. Flux densities and spectral indices for millisecond pulsars at 2.7 and 4.9 GHz. The 

difference in spectral indices derived for separate power laws below and above 1.4 GHz, 

Act = a<i.4GHz — o>i.4GHz is quoted in the last column. See text for details. 



PSR 


P 


-52695 


"54850 


spectral index 


diff. index 




(ms) 





mJy) 


(mJy) 




a 


Act 


J0437-4715 


5.757 






8.5 ± 2.5 a 


-1.17 ±0.06 


1.4 ±0.3 


J0621+1002 


28.853 


0.3 


±0.2 




-1.9 


±0.3 


6 


J0751+1807 


3.479 


0.8 


±0.3 




-1.2 


±0.2 


0.2 ± 0.8 C 


J1012+5307 


5.256 


1.2 


±0.2 


0.2 ± 0.1 


-2.0 


±0.1 


0.1 ±0.3 


J1022+1001 


16.453 


1.9 


±0.3 


0.53 ± 0.06 


-1.5 


±0.1 


0.0 ±0.5 


J1024-0719 


5.162 






0.4 ± 0.2 


-1.4 


±0.1 


d 


J1518+4904 


40.935 


0.8 


±0.2 


0.24 ± 0.07 


-1.5 


±0.2 


1.5 ± 0.8 e 


B1620-26 


11.076 


0.2 


±0.1 




-1.7 


±0.1 


1.5 ± 0.9 C 


B1640+2224 


3.163 


0.4 


±0.2 




-0.8 


±0.3 


.J 


J1643-1224 


4.621 


1.4 


±0.5 


0.4 ± 0.2 


-1.82 ±0.05 


0.5 ±0.5 


J1713+0747 


4.570 


2.7 


±0.9 


0.8 ± 0.2 


-1.2 


±0.2 


0.6 ±0.4 


J1744-1134 


4.075 


0.7 


±0.4 


0.19 ± 0.06 


-1.9 


±0.1 


0.4 ± 0.4 


B1855+09 


5.362 


2.3 


±0.7 


1.12 ±0.09 


-1.33 ±0.07 


0.2 ±0.3 


B1937+21 


1.558 


2.0 


±0.4 


1.0 ± 0.2 


-2.13 ±0.09 


0.1 ±0.1 


J2051-0827 


4.509 


0.3 


±0.2 




-1.5 


±0.1 


1.0 ± l.l c 


J2145-0750 


16.052 


2.2 


±0.5 


0.4 ± 0.1 


-1.88 ±0.08 


0.5 ±0.4 



a Measured by Manchester &: Johnston (1995) at 4.6 GHz. Uncertainty esti- 
mated to be 30%. 

b No flux densities at 1.4 and 4.9 GHz. 

c No flux densities at 4.9 GHz. Quoted value only based on 1.4 and 2.7 GHz 
measurements. 

d Formal fit suggests a flatter spectrum at high frequencies. 

c Possible turn-over around 400 MHz. 

f Only flux measurements at 0.4, 1.4 and 2.7 GHz available. 
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Table 4. Upper limits on flux densities for millisecond pulsars not detected at 2.7 and 4.9 GHz. 
Values are calculated for a 5a detection of that pulse feature which is most prominent at 1.41 

GHz (see Paper I). 



PSR 


P 


qupper 
'-'2695 


tyexp 
'-'2695 


Qupper 
'-'4850 


nexp 
'-'4850 




(ms) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


J0613-0200 


3.062 


0.46 


1.08 ± 0.60 


0.40 


0.52 ±0.31 


J0621+1002 


28.853 






0.32 


0.07 ±0.30 


J0751+1807 


3.479 






0.35 


0.40 ± 1.00 


B1257+12 


6.219 


0.38 


0.32 ± 1.00 






B1534+12 


37.904 


0.05 


0.09 ± 0.04 


0.04 


0.015 ± 0.007 


B1640+2224 


3.163 






0.44 


0.42 ± 1.10 


J2051-0827 


4.509 






0.44 


0.3 ±0.4 
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Table 5. Summary of polarization observations. We quote observing frequency, bandwidth, 
effective resolution and the degree of linear, n^, and total polarization, Ut- 





V 

(GHz) 


BW 

(MHz) 


Atcff 
(us) 


iiL ^ < I > 

(%) 




PSR 


n T = < vjw > 

(%) 


J1022+1001 


1.69 


28 


16.1 


32.7 ±3.1 


34.0 ± 3.1 




2.69 


40 


174 


16.6 ± 3.4 


18.0 ±3.1 




4.85 


40 


33.9 


3.7 ±3.8 


6.5 ±5.5 


J1713+0747 


2.69 


40 


6.4 


12.5 ± 1.3 


12.7 ± 1.3 




4.85 


80 


93.1 


15.3 ± 4.5 


16.7 ±6.0 


J2145-0750 


2.69 


80 


304 


5.9 ± 6.4 


12.3 ±5.1 




4.85 


200 


132 


3.1 ±2.2 


7.5 ±5.5 
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Table 6. Frequency dependences and values at 1 GHz of component separation (C), and profile 
width (W), respectively, and de-polarization index for MSPs with measured polarization at 5 GHz 

(see text for details). 



PSR 



Type 



Measure 
at 1 GHz (de£ 



7 



J0218+4232 


W a 


153 ±3 


-0.04 ± 0.03 




J0437-4715 


C b 


88.7 ±0.3 


0.00 ± 0.01 


-0.38 ± 0.05 


J0613-0200 


W c 


101 ±3 


0.00 ± 0.05 




J0621+1002 


C d 


104.0 ± 0.8 


0.00 ± 0.01 




J0751+1807 


Q d,e 


45 ± 1 


-0.02 ± 0.05 




J1012+5307 


W a 


47.3 ± 0.8 


0.00 ± 0.02 




J1022+1001 


cf 


12.3 ±0.7 


0.05 ± 0.08 


-0.56 ± 0.07 


J1024-0719 


C9 


27.3 ± 0.3 


0.0 ±0.1 




J1518+4904 


w h 


12.6 ± 0.5 


-0.03 ± 0.06 




B1620-26 


W c,d 


65 ±2 


0.08 ± 0.08 




J1640+2224 


w c 


49 ±2 


-0.02 ± 0.06 




J1643-1224 


w* 


31.3 ±0.8 


-0.27 ± 0.05 




J1713+0747 


w iJ ' 


13.5 ± 0.5 


-0.26 ± 0.05 


-0.34 ± 0.06 


J1744-1134 


w 


14.0 ± 0.6 


-0.2 ±0.1 




B1855+09 


C k 


23.8 ± 0.8 


-0.30 ± 0.03 




B1937+21 


w a 


13.5 ± 0.8 


0.0 ±0.1 


-0.76 ± 0.07 




c l 


172.6 ± 0.3 


0.00 ± 0.01 




J2051-0827 


w m 


26 ±2 


0.2 ±0.2 




J2145-0750 


w 


87.3 ± 0.6 


0.02 ± 0.01 


-0.5 ±0.2 



a Main pulse feature. 

b Using strongest outer components at —40° and +50° longitude. Additional data 
taken from Navarro et al. (1997). 

Additional data from Sallmen (1998) and Stairs et al. (1999). 

d Using strong outer components. 

"Additional data taken from Lundgren et al. (1995), Sallmen (1998) and Paper I. 

Using prominent leading and trailing component, see Kramer et al. (1999). 

s Using components remaining at 4.9 GHz. 

h Width referring to intensity of components remaining at 4.9 GHz. 

'Width measurements probably influenced by unresolved component at low frequen- 



j Additional data at 1.4 GHz taken from Sallmen (1998). 
k Using components as resolved at 4.9 GHz. 
'Using separation between peaks of main and interpulsc. 
""Additional data from Paper I and Stairs et al. (1999). 



32 



10* 

10 3 





10 




10 




10 




10" 1 




10" 2 


„ 




->, 




—3 


10 


E 






10 


->, 




-<— < 


10 


CO 




c 


10 


CD 




■D 


10 


X 




=3 


10 ' 


II 






10 




10 




10 




10 




10 




10 




10" 1 




10" 2 


„ 




->, 




—3 


10 


E 






10 


->, 




-4— ' 


10 


CO 




c 


10 


CD 




■D 


10 


X 




Z3 


10 


II 






10 




10 




10 




10 




10 




10 




10" 1 



icr 



J021 8+4232 i 



J0437-4715 i 



J061 3-0200 i 



J0621+1002 i 



-V 



J0751+1807 



J1 01 2+5307 



J1 022+1 001 



J1 024-071 9 



B1 257+1 2 



J1 51 8+4904 



B1 534+1 2 



B1 620-26 



J 1640+2224 



J1 643-1 224 



J1 71 3+0747 



J1 744_1 134 



B1 855+09 



B1 937+21 



J2051-0827 



J21 45-0750 



100 1000 100 1000 100 1000 100 1000 10000 

Frequency (MHz) Frequency (MHz) Frequency (MHz) Frequency (MHz) 



33 



J0218+4232 




50 -100 -50 50 

Longitude [deg] 



100 150 



34 



J0437-4715 



4600 MHz 




50 -100 -50 50 100 150 

Longitude [deg] 



35 



J0621+1002 




50 -100 -50 50 

Longitude [deg] 



100 150 



36 




50 -100 -50 50 100 150 

Longitude [deg] 



37 




50 -100 -50 50 100 150 

Longitude [deg] 



38 




50 -100 -50 50 

Longitude [deg] 



100 150 



39 



J1 51 8 + 4904 



4848 MHz 




50 -100 -50 50 100 150 

Longitude [deg] 



40 




50 -100 -50 50 

Longitude [deg] 



100 150 



41 




50 -100 -50 50 100 150 

Longitude [deg] 



42 




50 -100 -50 50 100 150 

Longitude [deg] 



43 



11855 + 09 



4850 MHz 




-150 -100 -50 50 

Longitude [deg] 



00 



50 



44 




50 -100 -50 50 100 150 

Longitude [deg] 



45 



J2145-0750 



4850 MHz 




50 -100 -50 50 100 150 

Longitude [deg] 



46 



1.0 




1 


1? 0.8 


_ J0751+1807 


M 2.7 GHz '_ 


1 0.6 


_ 


J \ - 


X 

M 0.4 




\ : 


Relative 

o o 
b k> 




v J M)h^ 




i 


' V W v ^ 

ii - 




150 200 250 
Longitude (deg) 



300 350 



47 



100 



80 



CC 




N 






60 


o 




Q. 




1— 




CC 




CD 




c_ 


40 


H 




O 




CD 




CD 




L_ 




D) 




CD 
Q 


20 







100 



• J0437-4715 
□ J1 022+1 001 
oJ1 71 3+0747 
^B1 937+21 
^ J21 45-0750 




1000 

Frequency (MHz) 



10000 



48 



CD 

■o 



c 
o 

CO 

1— 

CO 
Q. 
CD 
CO 



CO 
CD 
T3 



c 
o 

CO 

1— 

CO 
Q. 
CD 

co 



150 r 

100 7 

50 7 

o L 

100 - 

50 - 



u - 

50 r^ 

40 - 
30 - 
20 - 
10 - 

o L 

80 - 

60 - 

40 - 

20 - 



J021 8+4232 



J0621+1002 



J1 022+1 001 



=a=~Q~ 



B1 620-26 



CD 
T3 



^ 

50 r^ 

40 - 
30 - 
20 - 
10 - 

o L 



J1 71 3+0747 




J0437-4715 ~_ : 0613-0200 



J0751+1807 



J1 01 2+5307 



s O 



J1 51 8+4904 



J1 024-071 9 



J1 643-1 224 



- J 1640+2224 



J1 744_1 134 



Z 



i 



^r 



B1 855+09 



^=ir^ 



£ 100 - B1 937+21 

CO 

^_ 

CO 
Q. 

CD 50 - 
CO 



1 J2051-0827 



100 1000 

Frequency (MHz) 



J21 45-0750 



100 1000 100 1000 

Frequency (MHz) Frequency (MHz) 



49 



CD 
.Q 

E 

=3 



CD 
.Q 

E 

=3 



14 
12 
10 
8 
6 
4 
2 

18 

12 
10 
8 
6 
4 
2 




Normal pulsars 



Separation 
Width 



' I ' I ' I ' I ' I ' 


; msps 








- 


i i i i i i 








i 



-1.0 -0.8 -0.6 -0.4 -0.2 0.0 

Index 



0.2 



